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TORONTO HARBOUR NUMERICAL MODEL 
IMTRODUCTION 

As part of the overall Toronto Harbour study, a two-dimensional 
nuBierical model of currents and total dissolved solids (as 
iieasiired by conductivity) in the inner harbour has been 
developed. Such a model has its greatest value in its 
predictive capacity. The effect of water management decisions 
such as intake and outfall location, land filling and changes in 
harbour configurations and gaps can be simulated with a model 
of this nature. This is important as these projects involve 
large capital expenditures and considerable potential for water 
quality changes. 

Concentration contouring of physical-chemical parameters such 
as conductivity and dissolved oxygen can be developed for 
various discharge conditions, outflows, dispersion, 
Iteteorological conditions and changes in harbour geometry. 

For such predictions to be realistic, however j the model must 
be capable of considering phenomena which occur at very short 
time intervals (order of minutes). That such phenomena are 
important has been demonstrated in the case of Hamilton Harbour 
(MOE, 1974; Poulton and Palmer, 1976), where changes in 
conductivity of several hundred umho, and in dissolved oxygen 
of over 50% saturation, have been shown to occur within 
minutes. In the case of Toronto Harbour, Simpson (in Institute 
of Environmental Studies, 1974) states that short periodicities 
have little effect in enhancing water quality. With this 
conaideration in mind, recording chemistry and current meters 
have been operated for extended periods in the harbour and the 
real-time data thus obtained have been used directly in the 
numerical model. Internal points have also been instrumented 
and the results used for calibration and verification 



purposes. The present report describes the calibration of the 
model in terms of the physical constants used and the size of 
the time step required for stable solutions. Computerized 
contouring methods have also been applied for pictorial 
handling of spatial variations. 

DESCRIPTION OF MODEL 



With the exception of real-time input for the principal source 
and boundaries, the numerical model used is identical to that 
used in the Hamilton Harbour Study (MOE, 1974). The model is 
based on a two-dimensional estuarine water quality model 
initially developed and verified by J. Leendertse (1970, 1971) 
for Jamaica Bay, New York City, The depth-integrated equations 
for momentum, continuity and mass balance are coupled and 
solved by finite-difference techniques in a space-staggered 
grid. Although a three-dimensional numerical model of Toronto 
Harbour is presently under development, the initial work on 
Toroiito Harbour was conducted with a two-dimensional model. 
Such a model was felt to be adequate to describe the general 
circulation and water quality patterns in the inner harbour in 
the absence of thermal stratification (late fall); indeed, the 
presently available input data set is insufficient to calibrate 
and verify a three-dimensional model. However, the application 
of real-time data inputs to this model constitutes a real 
"first" as previous model studies generally used only mean 
daily loadings and hourly lake exchange values. Furthermore , 
the calibration studies involving values of the physical 
constants reported here will prove useful when the 
three-dimensional model currently under development becomes 
operational. 

The hydrodynamic (currents and water levels) portion of the 
model takes into account advection of momentum, water level 
variation, Coriolis force, wind stress, bottom friction and the 
effects on momentuii and volume of the addition or subtractioii 

of water at intakes and outfalls. 



The hydrostatic approximatioii is used and dispersion of 
momentim is ignored. The mass transport portion of the model 
takes into account advection, dispersion , and change of mass at 
intakes or outfalls. The equations and method of solution have 
been outlined (MOE, 1974) and will not be repeated here. 

Input Data 

A computational grid of space step length 152 m (500 ft) was 
prepared from a Canadian Hydrographic Service chart of Toronto 
Harbour, analogous to the Hamilton Harbour model. The grid 
used is shown in Figure 1. The east and west gaps were treated 
as open boundaries, while the Keating Channel (Don liver) was 
considered a source and the ship canal (Hearn generating 
station intake) a sink. To supply real-time data for input to 
the models the two gaps were monitored with digitally recording 
water chemistry and current meters; currents were recorded 
every 10 min, and conductivity every 20 min. The Keating 
Channel was monitored with a chemistry meter; as shipping 
consi derations precluded mooring a current meter in the 
channel, streamflow data from the Don liver gauge at Todmorden 
recorded at a 60 minute interval were used as input to the 
model. To check the validity of using flow data from one 
location and chemistry from another, the time series from the 
two locations were cross-correlated. Cross-correlation between 
streamflow and conductivity tended to indicate only a very weak 
relationship between the gauge and chemistry monitor. For 
example, two periods of increased flow (Figure 3) showed 
opposite corTelations with conductivity changes, likely related 
to runoff-caused dilution. Consequently, no time adjustment 
was used with these data. Locations of meter installtions are 
given in Figure 2; complete details of meter installations and 
operation are given in a separate report on meter operations 
(MOE, 1977). 

lo data were available for the Hearn ship channel on a 
real-time basis; indeed, average figures available varied over 
a tremendous range. Available flow values can be summarized as 
follows: 



S ource 

Simpson (Subcommittee Report) 

Hearn generating station, 

November 1975 

(H, Makuchj personal comm.) 

Lake Systems Unit 

(July 1976, unpublished) 



Flow (m s ) 



17-42; average 29 



144 



75-100 



With this large range of figures a¥ailable, it was decided that 

an acceptable figure would be one that allowed maintenance of a 
near zero change of average water level during the modelled 



period. Tests indicated such a figure to be 50 m s 



3 -1 



This figure may be sin,aller than the true figure as. storm water 
overflows were neglected in the first runs; the Toronto Central 
Waterfront Water Environment Study indicated several stom 

sewer overflows to the inner harbour with a total estiioiated 

3 -1 
flow of approximately 35 m /s (Subconmittee Report, p. 

13). Although these flows were not considered, the major 

portion of the storm water flow occurs in the Don River ^ which 

is considered in the model. On the other hand, as these flows 

represent a significant volume ^ the major flows should be 

monitored carefully in future surveys and included in the 

modelling efforts. 

Hourly values of wind data measured at the Toronto Island 
Airport were used as input to the hydrodynamic model. The wind 
stress coefficient relating wind speeds to surface stress was 
given the same value (0.0032) used in the Hamilton Harbour 
study (MOl, 1974). 



A modelling period was selected in which the three pairs of 
recording current and chemistry meters were all functioning 
well. The record also included two rainstorm events which 
produced significant loadings from the Don liver. The selected 



period was HoveBber 7-11, 1975, beginning at 1300 hours on 
November 7. The flow-conductivity relationship occurring in 
this period is shown in Figure 3. Hourly values for November 
7-8 used in the model are suamariEed in Table 1; current and 
chemistry meter data at 10-20 minute intervals though used in 
the ■odel are not included for the sake of brevity. In 
addition, the average wind values used are summarized in Table 

Selection of Model Parameters 

Initial attempts at operating the model in Toronto Harbour 
utilized nodel physical parameters taken from the Hamilton 
Harbour model study (MOE, 1974). However, as this was the 
first effort at operating a model with real-time data input at 
intervals of the order of minutes, and since large variations 
in currents and water chemistry are known to occur within 
minutes, it was felt tht a thorough investigation of model 
parameters was necessary for proper operation of the numerical 
model. This involved the checking of the hydrodynaiiic time 
step length, effect of wind stress on a moving surface, 
velocity techniques at impervious boundaries, depth and bottom 
roughness interrelations and the concentration time step 
length. 

Although the time step lengths of 2 or 2.5 min. as used in the 
initial Hamilton Harbour model runs produced nimierically stable 
results, Williams and Hinwood (1976) have indicated the need 
for a shorter time step to achieve replication of modelled 
results. Following work of Leendertse, they suggest the 
following criterion for the time step length t: 

1 _< At 'i/gE £2.5 

Ax 
where 

Ax = space step length 

g - acceleration due to gravity 

h - average or maximum water depth. 



Their own, woT'k suggested a value closer to 1, rather than 2.5, 
was required. For Toronto Harbour, for a water depth of 10 m, 
this criterion indicates that a time step length of the order 
of 15 to 37.5 s should be used. The nodel results should 
obviously not be a function of the time step length within this 
range. However, decreasing the time step length below 30 s 
produced rapidly increasing instabilities in which currents of 
over 100 ea a in either direction could occur within 
minutes of each other, clearly a physically unrealistic 
situation. In order to resolve this problem, the various 
forcing f'unctions of the niodel were examined in detail. 

As wind is the major forcing function for water movements , it 
was examined first. The wind vector should be a difference 
vector which considers surface water currents (J.B. Hinwoodj 
personal comm.). Model results after 18 hr for time steps of 
30, 15 and 7.5 s respectively (Figure 4) show some decrease in 
observed velocities but physically unrealistic values of up to 
72 cm s were still obtained. 

Examining the magnitude of the terms in the momentum equation 
at different modelled times revealed large variations in the 
advection terms near closed boundaries. This suggested 
alteration of the finite-difference approximations of these 
terms to replace velocities for points which represented land 
points, with adjacent water points. However, this procedure 
yielded no improvement in model results. 

On the other hand, success was achieved by increasing the water 
depth in shallow zones. As the wind stress term is inversely 
proportional to the water depth H, relatively high values of 
this term were found at points of very shallow depth, creating 
high velocities at these points which were gradually propagated 
to adjacent points. The model was therefore run with a minimum 
water depth of 4 m and time step lengths of 30, 15 and 7.5 s. 



A aample output of v velocities (velocity componetit parallel to 

Weaterti, Gap) as a function of time after 18 hr is given in 
Figure 5, the dotted lines again indicating zero velocity. 
Obviously, 30 aeconds is too long but 15 seconds was acceptable 
as it replicates the result at the shortest time step length. 
Next a time step length of 20 seconds was tested; the result 
(Figure 6) showed that this was also acceptable. Therefore, 20 
seconds was chosen as the longest acceptable tiiae step length 
in Toroiito Harboiir. 

However, the probl^ remained that the input data represented a 
cotisiderably deepened version of the harbour next to the 
islands along the southern and southwestern shores. In order 
to improve this situation, the miniJiiinD water depth was reduced 
to 2 n and the effect of bottom friction was next examined, 

Bott<m friction is ordinarily determined in the field by 
comparing a measured velocity profile to the wind stress 
results. Although some velocity profiles were measured in 
Haaiilton Harbour during 1976, the data were too variable to 

permit any meaningful calculation of Chezy's C or Manning's n . 

A literature survey of typical bottom friction equations was 

conducted; all figures were converted to Mannings n values in 

-1/3 
m s. Typical results are as follows: 



Location 

Mestemport Bay, 
Australia 

Lake Ontario 

Jamaica Bay, H,Y. 

San Francisco Bay 



Canal lined with coiicrete 0.018 

Irrigation canal, sand 
bottom 

Cobble-bottom channel 

River channel with 
even bottom 

Channel with irregmlar 0.074 
sides and bottom, with 
growth of weeds and grass 

Irregular channel with 0.119 
dense growths of willows 
and weeds 



Manning's n 


Reference 


0.052 


Williams & Hinwood 




(1976) 


0.05 


Simons (1971) 


0.034-0.038 


Leendertse and 




Grit ton (1971) 


0.130 


Leendertse and 




Liu (1975) 


0.018 


Chow (1959) 


0.030 


Ibid 


0.04,2 


Ibid 


0,052 


Ibid 



Ibid 



Ibid 



Without detailed field measurements, Manning's n must be 
estiMated by the fit of modelled data to observed data. 

Considering the range of values given above, and the bottom 

topography of Toronto Harbour, a grid consisting of four n 

values ranging from 0.040 (in the dredged part of the harbour) 

to 0.060 (in the shallow southeasterii part of the harbour) was 

established. This grid is shown in Figure 7. Although these 

values are arbitrary, they do reflect the known fact that 

shallow waters tend to have higher n values and provide at 

least an estinate for a parameter which is difficult to measure 

except in river systems . With this grid the maKimim velocity 

-1 



after 24 ho'urs was 12 cm s 
measured velocity of 11.3 cm s 



This is close to the maximim 
-1 



Effect of concent rat ion time step length 

A 24 hour model run on Toronto Harbour with a time step length 
of 20 a, requires 13 minutes of CPU time with 120K core on an 
IBM S/370-165 computer. The high cost of this computer time 
renders it difficult to model for an extended period of time, 
alter physical parameters (bottom roughEess, water depth, time 
step length, etc.) for model efficiency, or to test the various 
alternatives (harbour geOTetry, loading, etc.) that the model 
is designed for. Anticipated costs for Hamilton Harbour will 
be much higher due to the larger Bodelled area, unless the 
space step is increased to an unreasonably high value. 

The principal mechraism of mass distribution in the masff 
balance equations is advection. At the ma^cimimi velocities 
(15-25 cm s ) generally observed from 12 months of current 
meter data and a grid spacing of 152 m, one expects that a 
longer time step could be utilized in the model to save 
computer time. To this end, mass balance time step lengths of 
3, 5 and 10 times the selected hydrodynamic step length of 20s 
were tried. Siome sample results are given in Table 3 as time 
series of concentrations observed over a 2 hour interval at 
three points in the harbour. The mass balance time step of 200 
s was unacceptable as gradual cyclic changes in concentration 
were produced. The step lengths of 60 and 100 s both produced 
small random variations which were almost always below 0.5 mg/1 
except at the shallowest areas (Table 3b), where deviations of 
1*0-1.5 mg/1 sometimes occurred, especially at 100 s. 

It should be emphasized that the only justification in 
operating a mass balance time step longer than the hydrodynamic 
time step length is economics of computer time. The time step 
length chosen will represent a balance between accuracy and 
economics. Although operating at a 60 s mass balance time step 



length produces larger variations (+0.5 mg/1) which are larger 
than those produced when the hydrodynamic and mass balance time 
steps are equal (+0.1 mg/1), these variations are smaller than 
either measurem,ent precision (1 umho or 0.65 mg/1) or observed 
daily variations in the central harbour (daily standard 
deviations of conductivity in lovember 1975 were 2 to 5 umho, 
or a TD'S O'f 1.3 to 3 mg/1). If these observed variations are 
objectionable during cross correlation with field data, title 
series smoothing using binomial weights could be applied. 

The concept of increased mass balance time step has been 

successfully used by Wickramaratne et al (1976). They used a 

stability criterion based on advection as follows: 

max [ |u| + |v||J At < 1 

"Ax' ' 

Although they were able to use a mass balance time step of 15 
times their hydrodynamic time step, they were dealing with 
ocean currenta which are primarily tidally driven. In Toronto 
Harbour, prelioiinary hydrodynamic model results contained a 
principal periodicity of 8 to 10 minutes likely due to harbour 
oscillation. Cumulative effects from this oscillation are 
probably responsible for the deviations observed as the water 
quality tine step length is increased. 

MODEL RESULTS 

With the factors affecting model performance as discussed 
above, a ni^erical model run for 72 hours was performed , using 
input data of November 7-10, 1975. Pertinent input parameters 

used are sumnarized as follows: 

1. Average wind field according to Table 2. 

2. Wind stress coefficient 6 - 0.0032, 

3. Manning's w values according to Figure 7. 

4. Hydrodynamic time step - 20 sec. 

5. Mass balance time step - 60 sec. 

6. Minimum harbour water depth ~ 2 m. 

, . 2 -1 

7. Mass balance dispersion coefficient = 0'..'93 m""sec 
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As preliminary operations (Fi,gtires 5 tO' 6) indicated that 
significant periodicities O'f the order of minutes exist, 
■velocity (u atid v) , water level and concentration results were 
saved on tape every 2 min during the model run. 

TIME SERIES AMAL,YSIS Of MODBL PERIODICITIES 

A set of model u and v velocity data, at a 2 min interval from 

row 16, column 8 of the model grid (near current meter location 

129, Figure 2) was analyzed by standard time series techniques 

(Jenkins and Watts, 1968). The results are shown in Figures 8 

and 9. Most of the variance is concentrated in peaks of 12.8, 

8.4, 6.3s 4*7 and 4.3 min, which were generally significant at 

the 99% confidence level. Although some of the observed 

periods may represent aliasing of additional periodicities at 

periods of less than 4 min, these periodicities were in rough 

agreement with calculated harbour oscillation periods T 

n 

using the Merian equation for a symraetrical concave parabolic 
basin (Hutchinson, 1957, p. 304): 



Tn = '^^ 



(n[n+l]g Z^}h 



where : 

n = order of oscillation 

g - acceleration due to gravity 

z - maximiii. depth (10' m. was u,sed in calculation) 
£ = length (3200 n for longitudinal and 1700 m for 

t r ans ver ee os c i 1 1 a t i ons ) 

O'bserved and predicted periodicities are given in Table 4 
Model periods are separated on the basis that transverse 
oscillations have a stronger u velocity component and 
longitudinal oscillations a stronger v velocity component 
Differences between observed and predicted periodicities 
reflect the complexity of the actual basin. 



U 



MODEL PREDICTIONS OF CURRENTS AND WATER QUALITY 

Due to the presence of these osci Hat lows, half-hourly avexage 
values of currents and water quality were computed to give an 
overall trend. Predicted currents are shown in Figure lOA for 
conditions of predominantly west winds, and Figure lOB for 
conditions of predoniinantly east winds. The overall current 
picture in both cases shows the strongest currents following 
the wind direction along the shallow south shore, with a weaker 
return current occurring along the north shore. An overall 
counter-clockwise circulation is present under westerly winds ^ 
while an overall clockwise circulation occurs with easterly 
winds •, 

A computerized plotting representation of water quality 
contours predicted by the model is presented in Figures lla-lli 
for 6-hour intervals beginning from 24 hours of model time. An 
initial slow spread of the 225 mg/1 total dissolved solids 
(TDS) contour increases during the first portion of the period 
of influence of north-east to easterly winds (36 to 48 hours) 
with most of the harbour having TDS above 225 mg 1~ at 48 to 
54 hours. Subsequently, the advective effect of the return 
flow along the north shore combined with a strong inflow of 
lake water through the Western Gap causes a considerable 
retreat of the 225 mg/l area to cover only a small northeastern 
part of the harbour at 72 hours. The area of strong influence 
from the Don liver (defined by the 250 mg/l contour) remains in 
the northeastern corner of the harbour at all times. The 
influence of lake water entering the water at the Western Gap 
is shown by the 216 mg/l contour, ^lich reaches a maximum 
extent at 30 hours and retreats as the water with higher TDS 
advances across the harbour. This expansion and contraction of 
concentration contours is very significant as it demonstrates 
the non-steady condition of the harbour as shown by the model. 
The spatial variation of total dissolved solids is in 
qualitative -agreement with the average of a series of 1973 
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conductivity neasurements by the Great Lakes Division, 
Institute for Itivironniental Studies, University of Toronto 
(Fi,gure 12). These results show the lake water effect (less 
than 360 umho) in the northwest corner ^ the Don River effect 
(greater than 370 umho) in the northeast part of the harbour 
and relatively uniform values (363-370 umho) throughout the 
central and southern areas. 

VERIFICATION OF MODEL OUTPUT: CROSS CORRELATIOM 

In principle J model verification involves a large number of 
computer runs, each based on a different combination of input 
parameters such as surface and bottom stress, dispersion, etc. 
In order to present the capabilities of the present model, the 
variable parameters were s.-et as already described and a run of 
72 hr duration performed. Cross-correlation procedures similar 
tO' those used between current and water quality meter 
installations at the same point, or meter locations at nearby 
points (Poulton and Palmer, 1976), were adapted to a test of 
model validity. 

Three sets of cross-correlation were performed between the 
model and the meter data; these were: model u velocities with 
current meter north speeds; model v velocities with current 
m,eter east speeds; and model concentrations with conductivity.. 
Current meter directions were rotated sO' that north was 
parallel to the Y axis of the model grid (i.e. 55 east o^f 
true north) s and the left - handed model co-ordinate system was 
introduced into the cross correlation program so that all 
proper correlations would be positive. 

In addition to the model parameters already mentioned, four 
parameters can be varied which are pertinent to the 

croes-correlation processes. These are: 
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(a) Locatioti in model grid from which model data were taken. 

(b) Number of lags. 

(c) Humber of binomial sttoothing weights on model data. 

(d) Start time of meter data. 

With regard to model grid locations, eight locations were 
selected for cross- cor relation. These are shown in Figure 13. 
Locations A to F were selected for cross -cor relation with 
current meter data obtained at location 129; locations G and H 

were selected for cross-correlation with chemistry meter data 
at location 128. Preliminary attempts at current meter 
cross-correlation indicated a tendency toward some negative 
correlation;, as spatial plO'ts of mO'delled currents tended to 
support eastward flows close to the south shore and the 
eastward components predominated in the location 129 current 
meter data, locations D,, E and F were selected for current 
meter crO'SS-correlation. In addition, some averages of 
adjacent points were used. 

Most authors (foT example Shaatry, Fan and Erickson, 1972) 
reconmend that the miniinum number of lags be 10% of the series 
length; Fee (1969) suggests 10-15%. Mo'St examples quoted by 
Jenkins and Watts (1968) used a number of lags equal to about 8 
to iO'% of the series length. As the field time aeries 
corresponding to 72 hours of modelled data are quite short 
(chemistry J 144 data points; currents, 216 data points) it was 
decided to start slightly above lOl of the nuaber of data 
points and move downward if needed. For currents j 24 lags were 
selected and for chemistry, 20 lags were used. If the maximum 
cross -correlation occurs at a lag near zero, the coherence is 
recoiiputed with the crose-correlation function realigned to a 
lag of zero J as was done by Lee (1972) and Jenkins and Watts 
(1968). 

To reduce the effect of aliasing of short-term periodicities in 
the model data (Figure 8 and 9), the model currents (at a 2 min 
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conductivity measurements by the Great Lakes Division, 

Institute for Envircsnm.ental Studies,, University O'f TO'tonto 
(Figure 12). These results show the lake water effect (less 
than 360 umho) in the northwest corner, the Don liver effect 
(greater than 370 umho) in the northeast part of the harbour 
and relatively uniform values (363-370 umho) throughout the 
central and sO'Uthern areas. 

?ERIFICATIOM OF MODEL OUTPUT; CROSS CORIELATION 

In principle J model verification involves a large number of 
computer runs, each based on a different combination of input 
parameters such as surface and bottom stress, dispersion, etc. 
In order to present the capabilities of the present model, the 
variable parameters were set as already described and a run of 
72 hr duration performed. Cross-correlation procedures similar 
to those used between current and water quality meter 
installations at the same point, or meter lO'Cations at nearby 
points (Poulton and Palmer, 1976), were adapted to a test of 
model validity. 

Three sets of cross-correlation were performed between the 
model and the meter data; these were: model u velocities with 
current meter north speeds; model v velocities with current 
meter east speeds; and model concentrations with conductivity. 
Current meter directions were rotated so that north was 
parallel to the Y axis of the model grid (i.e. 55 east of 
true north), and the left - handed model co-ordinate system was 
introduced into^ the cross correlation program so that all 
proper correlations would be positive. 

In addition to the model parameters already mentioned, four 
parameters can be varied which are pertinent to the 

crO'SS-correlation prO'Cesses, These are: 
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Ca) Location in model grid froin which model data were taken. 

(b) lumber of lags. 

(c) lumber of binomial smoothing weights on model data. 

(d) Start time of m,ete.r data. 

With regard to model grid locations , eight locations were 

selected for cross-correlation. These are sho'Wn in Figure 13.. 
Locations A to F were selected for cross-correlation with 
current meter data obtained at location 129; locations G and H 

were selected for cro^ss-correlation with chemistry meter data 
at location 128. Preliminary attempts at current meter 
crO'SS'-correlation indicated a, tendency toward some negative 
correlation; as spatial plots of modelled currents tended to 
support eastward flows clO'Se to the south shore and the 
eastward components predominated in the location 129 current 
meter data^ locations D, E and F were selected for current 
meter cross-correlation. In addition, a'Ome averages of 
adjacent points were used. 

Most authors (for example Shastry, Fan and Erickson, 1972) 
recommend that the minimum, number of lags be 10i% O'f the series 
length; Fee (1969) suggests 10-15%. Most exOTaples quoted by 
Jenkins and Watts (1968) used a number of lags equal to about 8 
to 10% of the series length- As the field time series 
corresponding to 72 hours of modelled data are quite short 
(cheaiistry, 144 data, points; currents, 216 dat.a points) it was 
decided to start slightly above 10% of the number of data 
points and move downward if needed. For currents, 24 lags were 
selected and foT chemistry, 20 lags were used. If the maxim'um 
croes-eorrelation occurs at a lag near zero,, the coherence is 
recomputed with the cross-correlation function realigned to a 
lag of zero, as was done by Lee (1972) and Jenkins .a,nd Watts 
(1968). 

Xo reduce the effect of aliasing of shoTt-term periodicities in 
the model data (Figure 8 and 9), the model currents (at a 2 min 
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time interval) were snoothed with nine binoiiial weights and the 
chemistry with seven binoiiial weights prior to sub-sampling at 
the 20 minute (currents) and M minute (chemistry) frequenciee 
O'f field data. Variation O'f this param.eter within the range. 
above did not produce a large effect on the observed 
correlations , 

In some cases, maximim cross-correlations were initially 
obtained at lags other than zero. This is possible since the 
model interpolates wind data which are already 
.hourly-integrated values. In addition, the exact timing of 
current meter records (both input and verif iction.) is often 
uncertain. These digital records are sequenced by a 
fixed- interval timer and later converted to real time; hence, 
the times are approximate.. Evidence for the latter effect 
exists in the variable lags of the maximim cross-correlation 
results observed in the western gap during the model period 
(see report on Water Chemistry Meters in Toronto Harbour j this 
effect was allowed for in setting up input data to the model). 

The best results achieved in. cros.s-correlating model data with 
field data are summarized in Table 5. Plots of coherence 
sp€.ctra for u. and v velocities, .and concentrations appear in 
Figures 14-16j respectively. In each case, significant 
coherences were obtained for a portion of the accessible 
periods (greater th.an 1.5 hours for velocities and 2.0 hours 
for chemLstry) . 

Ideally,,, si,gnificant coherences should be obtain.e.d for the 
entire frequency domain for a perfect simulation. The 
crO'SS-correlation method is a very hard test for the model and 
the results indicate the importance of defining the model input 
constants such as bottom roughness and .surface drag, which are 
both temporal and spatial variables. In addition, the results 
emphasize the need for metered data (both input and 
verification) of high quality., Observed data at location, 129 
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indicated zero speeds abo^ot 30'% of the time and sudden changes 
from over 10 cm s ' to zero. In addition, many of the 
cofflpass readings in the eastern gap indicated a flow 
perpendicular to channel which is physically unrealistic. To 
maintain constant water level in the harbour during modelling 
it was necesaary to assime that this condition represented an 
outflow of water from the harbour.. 'This has undoubtedly 
affected the performance of the mod,el with the present data. 
However, it must be remembered that the maintenance of 4 or 5 
pairs of digitally recording current and water quality meters 
in simultaneO'US accurate operation is in itself a difficult 
task, one that must be accomplished properly in order to 
operate models of this type. 

If the model is working properlyj the spectral analysis of 

model data and measured data carried out separately should show 
similar periodicities. Even with smoothing as carried out in 
this analysis J the model spectra are heavily contaminated with 

aliasing of high frequency peaks (see Figures 8 and 9) when 
sampled at 20-30 minute intervals. For realistic comparison of 
model data with field data, a continous recording current meter 

is desirable. Use of a hot film anemometer to derive 
high-fre'quency kinetic energy spectra, as has been done in lake 
Ontario ( Palmer, 1973) is desirable. Given data of this 
nature,, adjustment of model physical parameters and 
verification of the model should become less of an onerous 
task. 
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TABLE li Model input data, November 7-8, 1975. 



XD 







East 


Gap 


West 


Gap 


Keating 


Channel 


Wind 








Current 


Cond . 


Current 


Cond . 


Flow 


Cond. 


Speed 


Direction 


Date 


Time 
13.00 


cm/sec 
1.40 


umho/cm 
333 


cm/ sec 
1.37 


umho/cm 
320 


mS/eec 
1.59 


umho/cm 
900 


m/sec 
3.1 


(from) 


Nov 7 


ESE 




14.00 


0.68 


333 


2.47 


320 


1.59 


900 


4.0 


S 




15.00 


0.67 


334 


-0.37 


320 


1.61 


900 


3.1 


ISE 




16.00 


4.03 


333 


-3.97 


320 


1.68 


900 


2.7 


SSE 




17.00 


2,90 


335 


4.26 


319 


1.75 


828 


2.7 


E 




18.00 


1.11 


334 


7.90 


319 


1.86 


746 


0.4 


SW 




19.00 


1.09 


337 


6.63 


320 


1.99 


718 


5.4 


SSW 




20.00 


0.63 


336 


6.07 


320 


2.30 


616 


6.7 


SSW 




21.00 


1.38 


335 


8.20 


320 


3.71 


814 


3.6 


WSW 




22.00 


4.15 


338 


8.43 


320 


5.07 


710 


1.8 


WSW 




23.00 


0.41 


338 


5.64 


321 


4.96 


850 


2,7 


SSW 


Nov 8 


0.00 


1.28 


336 


-2.88 


321 


4.31 


892 


2.2 


WSW 




1.00 


0.83 


335 


-11,28 


321 


3.80 


976 


1.3 


SSW 




2.00 


1.31 


334 


-5.62 


321 


3.40 


974 


2.2 


SSW 




3.00 


0.07 


334 


4.08 


320 


3.03 


9 62 


4.5 


SSW 




4.00 


0.49 


335 


-3,70 


320 


2.67 


976 


2,2 


SW 




5.00 


0.39 


335 


-5,05 


320 


2.45 


922 


2.2 


SSW 




6.00 


0.09 


335 


2.76 


321 


2.36 


928 


1.8 


S 




7.00 


0.07 


333 


5.80 


322 


2.27 


9 34 


5.4 


s 




8.00 


0.07 


333 


7.98 


320 


2.18 


924 


3.1 


SW 




9.00 


0.35 


333 


5.95 


320 


2.08 


582 


4.0 


SW 




10.00 


1.51 


333 


2.65 


321 


2.04 


926 


5,4 


SW 




11.00 


0,45 


333 


2.61 


319 


2.01 


818 


6.3 


SW 




12.00 


0,55 


333 


2.43 


320 


1.97 


838 


6.3 


SW 




13.00 


1.18 


331 


-4,72 


319 


1,95 


844 


6.7 


SW 




14.00 


1.82 


332 


5.39 


319 


1.97 


690 


7,2 


w 




15.00 


2.37 


336 


3.78 


320 


2.05 


814 


8,5 


w 




16.00 


-2.83 


336 


-2.99 


318 


2.14 


702 


8.5 


w 




17.00 


-2.15 


332 


3.22 


317 


2.21 


794 


7.6 


WNW 




18.00 


-3.23 


332 


4,49 


318 


2.15 


804 


7.6 


WMW 




19.00 


-2.76 


334 


7.49 


328 


2.09 


814 


6.7 


WNW 



Note: Negative current indicates flow out of harbour. 



TABLE 2 : Average wind speeds encountered during 7 2 -hour 
period modelled 



Date 

Nov 7 
Nov. 7-8 
Nov. 8 
Nov. 8 

Nov, 8=9 
Nov. 9 
Nov. 9-10 
Nov. 10 



Time 



13.00 

19.00 

7.00 

14.00 

20.00 

2.00 

9.00 

9.00 



16.00 
6.00 

13,00 

19.00 
1.00 
8.00 
8.00 

13.00 



Av, direction 
(from) 


Av . Speed 
in./ sea) 


SSE 


3.2 


SW 


3 


sw 


5.3 


W 


7 . 7 


w 


t 


NE 


2.4 


m 


6.3 


i 


8.1 
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TABLE 3 : Model concentrations for 2 hours at various concentration 
time step lengths. 

(a) Central deep location (about 10 m deep) 

Concentration time step/hydrodynamic time step 



lOX 


5X 


m 


IX 


0.2099 


0.2098 


0.2090 


0.2096 


0,2 097 


0.2094 


0.2096 


0,2096 


0.2101 


0.2099 


0.2098 


0.2096 


0.2102 


0.2099 


0.2098 


0.2096 


0.2102 


0.2098 


0.2098 


0.2096 


0.2104 


0.2098 


0.2097 


0.2096 


0.2105 


0.2098 


0.2097 


0.2096 


0.2102 


0.2097 


0.2097 


0.2096 


0.2105 


0.2095 


0.2095 


0.2096 


0.2105 


0.2097 


0.2098 


0.2096 


0.2103 


0.2096 


0.2096 


0.2096 


0.2100 


0.2097 


0.2097 


0.2096 


0.2101 


0.2097 


0.2097 


0.2096 


0.2099 


0.2097 


0.2097 


0.2096 


0.2099 


0.2100 


0.2099 


0.2096 


0.2095 


0.2095 


0.2095 


0.2096 


0.2099 


0.2098 


0.2097 


0.2096 


0.2094 


0.2096 


0.2095 


0.2096 


0.2094 


0.2097 


0.2097 


0.2095 


0.2094 


0,2 095 


0.0295 


0.2095 


0.2093 


0.2097 


0.2097 


0.2095 


0.2091 


0.2098 


0.2098 


0.2 095 


0.2089 


0.2 098 


0.209 8 


0.2095 


0.2089 


0,2098 


0.2 096 


0,2095 


0.2088 


0.2099 


0.2098 


0.2095 


0.2007 


0.2098 


0.2097 


0.2095 


0.2088 


0.2097 


0.2097 


0,2 095 


0.2089 


0,2096 


0.2095 


0.2095 


0.2088 


0.2096 


0.2095 


0.2095 


0.2091 


0.2098 


0.2098 


0,2095 


0.2089 


0.2095 


0.2096 


0.2 095 


0.2092 


0.2097 


0.2 097 


0.2095 


0.2089 


0.2096 


0.2 096 


0.2095 


0.2091 


0.2097 


0.2096 


0.2095 


0.2090 


0.2095 


0.2095 


0.2095 



3 
Mote: All concentrations are in kg/m (mg/lTlOOO) , and printed 

at an interval of 200 sec. Hydrodynamic time step is 

20 sec. 
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TABLE 3 : co^ntinued 



Model concentrations for 2 hours at various concentration tiiie 
step lengths. 

(b) Southern shallow location (about 2 m deep) 

Concentration time step/hydrodynamic time step 

lOX »X ^ 3X IX 



0. 


2093 


0. 


2097 


0, 


2091 


0. 


2087 


0. 


2093 


0. 


2096 


0. 


2082 


0, 


2099 


0. 


2088 


0. 


2086 


0. 


2086 


0. 


2066 


0. 


2087 


0. 


2051 


0. 


2066 


0.. 


2063 


0. 


2059 


0. 


2060 


0. 


2064 





.2063 





.2079 





.2062 





.2086 





.2082 





.2089 





.2108 





.2101 





.2128 





.2125 





.2144 





.2143 





.2155 





.2156 





.2156 





.2154 





.2171 



0. 


2099 


0. 


2098 


0. 


2098 


0. 


2101 


0. 


2097 


0. 


2109 


0. 


2091 


0, 


2018 


0. 


2098 


0. 


2100 


0. 


2105 


0, 


2 095 


0. 


2112 


0. 


2 097 


0. 


2109 


0, 


2106 


0. 


2104 


0, 


2105 


0. 


2107 





2097 





.2112 





.2093 





.2111 





.2099 





.2102 





.2107 





.2092 





.2106 





.2095 





.2097 





.2097 





.2096 





.2095 





.2095 





.2086 





.2103 



0.2102 
0.2096 
0.2099 
O.?103 
0.209 3 
0.2107 
0.2090 
0.2105 
0.209 7 
0.2098 
0.2101 
0.2093 
0.2103 
0.2101 
0.2100 
0.2102 
0.2097 
0,2098 
0.2101 
0.2092 
0.2106 
0.2089 
0.2104 

0,2096 
0.2099 
0.2099 
0.2093 
0.2100 
0.2099 
0.2096 
0.2101 
0.2094 
0.2097 
0.2099 
0.2089 
0.2104 



0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0-2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0-2099 
0.2099 
0.2099 
0.2099 
0.2099 
0.2099 
0,2099 
0.2099 
0.2099 
2099 
2099 
209 8 
2098 
0.2098 
0.2098 



Note : 



All concentrations are in kg/m (mg/lTlOOO) , and 

printed at an interval of 20 sec. Hydrodynamic 
time step length is 20 sec. 
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TABLE 3 : continued 

Model concentratione for 2 hours at various concentration time 
step lengths. 

(c) Southwestern location (about 8 m deep but close to shore) 

Concentration time step/hydrodynamic time step 

lO'X 5X m IX 

0.2099 0.2100 0.209 7 

0.2092 0.2094 0.2097 

0.2100 0.2100 0.2097 

0.2094 0.2096 0.2097 

0.2094 0.2095 0.2097 

0.2095 0,2097 0.2097 

0.2 093 0.2096 0.209 7 

0.209 8 0.2098 0.209 7 

0.2093 0.2096 0.2097 

0.2094 0.2095 0,2097 

0.2098 0.2099 0.2097 

0.2092 0.2095 0.2096 

0.2094 0.2096 0.2096 

0.2091 0.2096 0.2096 

0.2092 0.2096 0.2096 

0.2094 0.2098 0.2096 

0.2087 0.2094 0.2096 

0.2095 0.2099 0.2096 

0.2092 0.2097 0.2096 

0.2092 0.2095 0.2096 

0.2093 0.2098 0.2096 

0.2090 0.2095 0,2096 

0.2095 0.2096 0.2096 

0.2089 0.2095 0.2096 

0.2089 0.2094 0.2096 

0.2094 0.2099 0.2096 

0,2090 0.2095 0.2096 

0.2093 0.2098 0.2096 

0.2092 0.2097 0.2096 

0.2093 0.2096 0.2096 

0.2095 0.2099 0.2096 

0.2088 0.2093 0.2096 

0-2094 0.2097 0.2096 

0.2090 0.2094 0.2096 

0.2092 0.2093 0.2095 

0.2092 0.2095 0.2095 

3 
Note: All cO'ncentrations are in kg/m, (m,g/lTlOOO) ,, and 

printed at an interval of 200 sec. Hydrodynamic 

time step length is 20 sec. 






.2097 





.2090 





2096 





.2091 


0. 


2094 





2093 


0. 


2092 


0. 


2098 


0. 


2089 


0. 


2091 


0. 


2097 


0. 


2093 


0, 


2097 


0. 


2092 


0, 


2098 


0. 


2100 


0, 


2093 


0. 


2104 


0. 


2099 


0. 


2101 


0. 


2102 


0. 


2100 


0, 


2108 


0. 


2102 


0. 


2102 


0. 


2105 


0. 


2101 


0. 


2105 


0. 


2101 


0. 


2103 


0. 


2108 


0. 


2100 


0. 


2106 


0, 


2104 


0. 


2104 


0. 


2100 



'm 



TABLE 4: Toronto Harbour Oscillation Periods as Deteirmined 





by Numerical 


Model 


and Merian 


Equation 




(a) 


Longitudinal 












Model Period 


(min) 


Merian Eqn. 


Period 


(min) 


1, 


12.8 






12.0 




1 


8 . 4 






6.9 




1 


4.7 






4.9 




(b) 


Transverse 










a 


Model Period 


(min ) 


Merian Eqn, 


Period 


(min) 


1 


6.3 






6.4 




1 


4.3 






3.7 





Note: Model periods are obtained from spectral analysis of 
model u and v velocities. 

Merian equation periods are calculated from the equation 
for a symmetrical concave parabolic basin (Hutchinson, 
19 57, p. 304) . 
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TABLE 5 : Significant Coherences between Model Data and Field 
Data, Toronto Harbour 



Data Model Location Meter Location Significant Periods, hr 

U velocity Average of A, C 121 1.8-5.3 

V velocity 'W lii 1.1-2.0 

5.3-16.0 

Concentration G 12i 1.4-1.7 

(Dissolved 2,2-2.5 

solids) 6.7 

20.0 

Mote: Components of current meter velocity parallel to x and 

y axes of model grid (i.e. east and west gaps) were used 
in cross-correlation with model U and V velocities, 
respectively. 
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FIGURE I -.TORONTO HARBOUR NUMERICAL MODEL GRID 
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FIGURE 2 : LOCATION OF BOUNDARY AND INTERIOR MONITORING POINTS FOR TORONTO HARBOUR NU 



MERIGAL MODEL 
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riGURE 3^ RELATIONSHIP BETWEEN DON RIVER STREAMFLOW AND KEATING CHANNEL CONDUCTIVITY NOV. 7-11,1975 
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FIGURE 4 -ORIGINAL PREDICTION OF NUMERICAL MODEL CURRENTS AT TIME STEPS OF 30,15 AND 7-5 SECONDS 




FIGURE 5 '= PREDICT lOiNiS O'F NUMERICAL MOD€L CURRENTS WITH REVISED WATER 

DEPTHS AT TIME STEPS OF 30„ 15 AND T-5 SECONDS 
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FIGURE 6 : PREDICTIONS OF NUMER ICAL MODEL CURRENTS WITH REVISED WATER 
DEPTHS AT TIME STEPS OF 20 AND' 15 SECONDS 
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FIGURE 7 ' manning's N VALUES USED IN TORONTO HARiOUiR ■HODEL 
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FIGURES: AUTOSPECTRA OF UNSMOOTHED MODEL U VELOCITY AT EAST METER LOCATION (129) 
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FIGURE 9 ■ AUTOSPECTRA OF UNSMOOTHED MODEL V VELOCITY AT EAST METER LOCATIOH (129) 
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FIGURE 10 a ^ NUMERAL MODEL CURRENTS AFTER 27 HOURS (WIND PRIMARILY FROM 
WEST, AVERAGE SPEED 8 M/S) 
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FIGURE 10' b : NUMIERICAL MIODEL CURRENTS AFTER 80 HOURS (WIND' PRIMARILY 
FROMi EAST FOR PREVIOUS 24 HOURS, AVER SPEED 6 M!/S) 
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i''xgure 



11a 



iJuKierical rLOdel PredictiO'ns of 'VO'tal Ddsoolved Solids 
After 24 Hours of Modelled '.Jirie 



Indicated figures are values O'f total dissolved 
solids (mg/l) for each contour line. 
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Figure lib 

Mumerical Ilodel Predictions of ''JO'tal Dissolved Golids 

After 30 Hours of riodelled '^^iue 



Indicated figures are values of total dissolved 
solids (iEg/1) for each contour liiie. 
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Figure Uc 

lumerical Model Predictions of Total Dissolved Solids 
After 36 Hours of Modelled Time 



Indicated figures are values of total dissolved 
solids (mg/1) for each contour line. 
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Figure lid 

Numerical ilodel Predictions of Total Dissolved Solids 

After 42 Hours of node lied Time 



Indicated figures are values of total dissolved 
solids (mg/l) for each contour line. 
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Figure lie 

Numerical nodel Predictions of Total Dissolved Dolids 
After 48 Hours of Ilodelled Tine 



Indicated figures are values of tO'tal dissolved 
solids (mg/1) for each contour line. 
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Figure llf 

numerical llodel Predictions of Total Dissolved Solids 
After 54 Hoors of Ilodelled Time 



Indicated figures are values of total dissolved 

solids (mg/l) for each contour line. 
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Figure llg 

Numerical Ilodel Predictions of Total Dissolved Solids 
After 60 Hours of Ilodelled Time 



Indicated figures are values of tO'tal dissolved 

solids (mg/1) for each contour line. 
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Figure llh 

Numerical Ilodel Predictions of Total Dissolved Solids 

After 66 Hours ox Ilodelled Time 



mm 



Indicated figures are values of total dissol¥ed 
solids (mg/l) for each contour line. 



tVillLHiJtkCUiLLLltlLULiuU.llttiM..>.. 



XttltltlUittlttyi 



^tfeiittU 



'I'iiLIL.Ik.LLil.M 
TUliilLi.tll.ltlLIH.f4tititiiife.L'L!L 



t kUl till Li!i LILU L^ LLIL i Lti 

|UlilltUU(.|iUL<IL.II.UIL.LL%iLlt[L.l.).ilL 



4ILUi 



"<L!tUilliLiL(.LlllLILtL<U i.i lUlLilLlLCLLL 

'UUlliLLiLtill,ltiLtUtU.aULUilLILtU 

ititLiLUtLaiuikuitii!iLLVll.ibtUyli 

.■IHLtilkULflttL'ILtL 

XLtUiiLLILlLiUt 

.i!Ui'!LttlLitiUl.'ULllLiULlllLl.illkiLlt 

lHl^.^|■|i!iu.^lltt^uiJitiw^i;tluulLl.l.U 



.ll=|iLLUtitLiM.iiUi 



™».Ll,lH,l.ilt 
ILUILliUl-l.' 



ittUllLlLktiLLU 



lLLIUUILHIlLL!k,L 



ILLIUiUUiLSiUULtU 
(LiitLUiLlUULLULklLi 

>UkLUj;ii|i,LUlllLn;Lil.L 

U,tl.t.ilLttl).ili:l"i!4,tl.lttL 
ILUt-LLLaLILILitUi. 

.lUklitliiL 



ii;4LLiUULLi. 

t».tl.l.l.l.iLtLL«.it.iyi.il:t 
H,LULILilb.iU),ui.iii 



*■'* -i L- K" **" 



2 26 



2 2^.,...:: 



216 ::: 




I 
I 



t 
I 
I 



216 



I aicij,i,4iiuiiLiiL:aCi,L> 



..ILiiliiLJ 
.4i.iLU.it I 



.LiLkLLI 
.Ltttt'M 



tiLi!.l.iitll.Lli.lLLI.iLLLIL9.l I 



,.lLtLll,!tlLI 
'LLtULiLLl 



,iLl,iLi,L'iLLLi,ii.iULIk1i.^lL^tU' 



X'U il.Lll.lL4l,0LLUaL!LILU.I 



«LU.LUl 



LLlLL'ttLUllUi 



4f''L«,Ul.i 

.'|.4Lit.ltiil.iCll.iLLL.Ui 



.ULI^LLiLLLIiiLLtiLU 

L'LL LU Hit til L LI LLiUia. I. U L 
■AtWLtl.iU4.U.l.LtlUl.U« LLiLL 

.iLULLLlL.LLILLL.|kltLiti.l 

■ " titl.ilL.ll.iLLL.ILIl 

UtLiULHLiI, 

.LtliiiLU. 
■ LitLUlti 



4LUILIL 



.titLLIHtlJlLtittLtiittttl-iiULI 
,tlL,iLLLiitt*'LtLUtlUl.LLiLtLlt.lL;| 

.llill.l,tLLtLlLii,l,LUUtULtULi| 

ttiiUttl 
.LLilLilUUUIkU.LL.ilUU^iUULiuLiikllUta 



.iLiLLLl 

ILiitiiLl 

HLi)LLit.Lfc| 



.l.iLtLiil.ltLl.tlL{ 



i(.LiULILLLilLLilLL'LiLl.ULl.Li 
LtlULULibU^iiiitULLLt 
4tiUML|iLi,t Li 

tl,ll.lljLtiLLlnl 

iLRiiLk<.i 



L4i I tLfeiife.i.|,tt ifciL L lit Li 
iiSiULlLilLlli».iiill.iLLLL ,„=. 

i.LLiuutitLktkLakLii.i,iuiuLiki 



.liiLLii 



.iLLHiLi 



<L.lLiiU 



XtLll.tL.iL(,l.LI.' 



,lLil..LltLLiLlLlLiLLILi 



.itl.LltU'I.UUILiiLLiitHililLl,iLLiUILi 
jUtLlH.!L.ILlil.LILHi.itH.ILliLUUtil 



U^LLlLLiLLl,' 

iUbiLLLULi 

.H,ttiH.(LLLLilil.t! 

tltLiliuiLLLi 



LiLiLi^iiiLLii 

iiLLLiLLiaLiLLL,ILLU.tLl,ULin,LHii.itliillLii 
iL LIL4 LyibU Lit LUiLLiLLi il Lit L LiL L.^IL Li IL kUl 
iLUiLI.Li^Li.LLill.L:ILULilL'D]LlL'LlU4LL<Ui..Lil 

jLUiiLU,H,iiLiLiiLi,iU|.ikltiH.l,iLLlLii 

.iLUIH,tiULB.iLLIkitLiilLiL.tiiLLtiLttil 

.ILiLiLi.l,iLLIliLLiLll.iLLLi«.HiLLII.Ul'IH.I 



I 
I 
I 
I 



.i'1,tfeiii|tllULi!Ltt.i 
JitttlUilfUiH 



44 



I 
I 



Figure Hi 

Huifierical Ilodel Predictions of Total Dissolved Solids 

After 72 Hours of node lied Time 



Indicated figures are values of total dissolved 
solids (mg/1) for each contour line. 
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FIGURE 12 : AVERAGE TORONTO HARBOUR CONDUCTIVITIES (Mmho) MEASURED IN 1973 




FIGURE 13 • LOCATION OF MODEL DATA USED FOR CROSS CORRELATION 
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FIGURE 14- CROSS CORRELATION OF MODEL U VELOCITY (AVERAGE OF POINTS A AND C) WITH COMPONENT OF 
CURRENT METER DATA PARALLEL TO EAST GAP, LOCATION 129. 
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FIGURE IS: CROSS CO'RRELATIO'N OF MODEL V VELOCITY AT LOCATION ID WITH COMIPO'NENT OF CURRENT METER DATA 
PARALLEL TO WEST GAP, LOCATION 129. 
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APPENDIX 

TWO-DIMENSIONAL NUMERICAL MODEL COMPUTER PROGRAM 

The computer program as used In this study is capable of 
solving the coupled equations of momentum, continuity, and mass 
balance as described by Leendertse and Gritton (1971). 
Although so far unused, the program Is capable of modelling up 
to five coiiponents , and provision is made for first-order decay 
or interaction rate constants for these components. The user 
is capable of selecting previously created real-time input 
files to describe flows and concentrations at sources and open 
boundaries, or to enter averaged data at model run time for 
these parameters. An update feature permits restarting of the 
program from any point in modelled time of a previous run. 

The purpose of this Appendix is to describe the required input 
data files and provide a sample of model output, for the 
assistance of new users unfamiliar with the program. 
Descriptions of two auxilliary programs used to produce 
real-time input files for the model are also included. Copies 
of the programs (listings, cards and/or magnetic tape) are 
available on request from Dr. D. Poulton, Water Resources 
Branch . 

The program is written in FORTRAN IV compatible with IBM 360 
and 370 installations, and is composed of the following 
routines: 



MAIN: establishes common blocks required for data storage 
for any geographical location. 

HARBOR: reads in all physical data, controls calling sequence 
to all other subroutines and program output. 

INPUT; enters real-time data (winds, source and boundary 

flows and concentrations) to the appropriate files. 



A-1 



FRITE: produces outpot maps of selected variables 

(velocities, water levels, and concentrations) at 
specified intervals. 

IDB; calculates double-bar water depth means. 

CHEZYi calcolatee Chezy coefficients. 

UVIL: calcolates U velocities and water levels at half-time 

steps . 

VVIL: calculates ¥ velO'Clties and water levels at full-time 

steps. 

PCO'Hi calculates cO'nstitu,ent concentrations at half-time 
steps . 

QCO'N: calculates constituent concentr at io^os at full-time 
steps . 

Tape or disk files used by the model are given in Table A-1.. 
Pile 10 is model output, whereas the other files are input. 
File 20 is presently created internally, while files 25 to 46 
■ay be either created internally or copied from external tapes 
before model operations. All files are unformatted? files 10 
and 99 noraally reside on magnetic tape and the other files on 
disks, though this may be altered through the Job Control 
Language . 

A complete descriptiO'n of the input data cards required for 
running the model prograiii is given in Table A-2. Mote that 

there are 32 types of data cards? more than one of certain 
cards is required, and if real-time sO'Urce and boundary data 
are specified, certain types o^f data cards are omitted 
entirely. The fO'rm.ats o^n many of the cards co'ntain blank 
fields which are useful for writing identifying information on 
the cards. This fact may be clarified by reference to Figure 
A-1, which is a eo^mplete listinf of the data deck required to 
run the numerical model as used in this report starting at 24 
hours and ending at 72 hours.. Note that NSECW'L, NSECF, and 
NSECC are and consequently cards type 20,, 21^ 22, 24,. and 26 
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do not appear. Card 29 must be present even though its value 
is not used. Although not read by the model, columtis 10 to 16 
of the wind cards contain the year, month, and day of wind data 
followed by •" 1" or "2" for first or second half of day. 

Two additional programs have been written in this study to 
transform current and water quality meter recording data to the 
proper format accepted by the model, one for open boundaries 
and one for sources and sinks. 

Input data used with the open boundary program are given in 

Table A-3. The first two data cards define the modelled period 
and allow for separate starting times for the current and water 
quality meter data records, in case there is an error in the 
time base for either set of recorded data (detected by a lag 
time to maximum cross-correlation; see report on Toronto 
Harbour Water Chemistry Meters). The number of spinup steps is 
that used in the model itself (Table A=2, Card 3) divided by 
the interval between successive records for real time data 
(NSTEPS, Table A-2, Card 6). The model time length, MINT used 
here is the hydrodynamic time step length multiplied by this 
same factor NSTEPS; this factor is used to reduce computer time 
for data input and is normally set so that MINT' is the highest 
common factor in MINC and MINQ. The program decides the sign 
of the current meter speeds according to the compass value 
compared to CHANEL; if the compass reading is ambiguous (more 
than 45° from CHANEL or CHANlL+180), the total dissolved 
solids (conductivity x 0'.65 from the chemistry meter) is used 
to determine the sign.. Assuming that harbour water has a 
higher total dissolved solids (TDS) then lake water, values of 
TD'S above a maximum value TDSMAX produce a negative (out of 
harbour) current and vice versa. The total dissolved solids 
data are also checked against the currents; if the TDS is above 
TD'SMAX for a positive current, TDS is set to 0.9*TDSMAX. The 
speeds and dissolved solids are interpolated to the model data 
file time step length if necessary, and output to tape files 
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corapatible with model files 35, 36, 45 ana 46. In addition, 

printed output is produced, one line for each outpat record,- a 
star is printed to indicate an uncertain compass value, with 
the sign of the speed set by TDS , or the TDS being limited by 
positive speed. It should be noted that the present program is 
limited to model variables NCHANR, NCHANC, and LMAX all ^ 1 and 
must be modified for other conditions. 

The source-sink program is similar to the open boundary 
programi, except that it allows for up to 16 sources or sinks , 
some of which may have only average (constant) or slowly 
varying (e.g. dally) flow and concentration data. Input data 
formats are given in Table A-4, The function of spinup steps 
and time step length is identical to the boundary pro^gram, i the 
sign convention for the current meters is similar, except that 
there is no cheek according to the dissolved solids value » 
Current m,eter speeds are multiplied by the croes-seet ion area 
of the channel or slip at the meter location to get a flow in 
M 's , and all flows are divided by the area of the 
grid-size to give a flow speed (m s~ ) into the grid square 
In which the source is located. The data are inter po'lated to 
the ii;Odel file time step length as needed, and output to tape 
files compatible with model files 2 5 and 30. At present this 
program is also limited to LMAX ^ 1. 

The regular version of these programs assumes that currents and 
water quality are on tape files of standard format as used in 
the Lake Systems Unit, Fat currents, this means formatted 
blocks of 24 0' words, each written in format 24 0F10'., 2, 
containing lines of five words representing date, time, 
temperature,, compass (from,) and speed (cm, s~^). For 
chem,istryy this means unformatted blocks o^f 324 words each 
containing lines of nine words with date and time being the 
first two words and conductivity the last word. The other 
words of each line are unused,. 



For Toronto Harbour only, the source program was modified to 
replace current meter input with digitized Don River flows 
sopplied by the Water Survey of Canada on data cards with 
foriiat 16X,8F8.2. 

Model Outpu t 

The nunerical model produces two types of printed output and 
one type of magnetic tape output (File 10, Table A-1), A 
sample of output produced at each time step by routine HARBOR 
is shown in Figure A-2. Output consists of the time step 
number, water volume above datum in m , wind speed (m s ) 
and direction (with respect to model y-axls), row and column 
channel currents in m s , first and last source flows and 
concentrations in m s and kg m respectively, and 
output velocities (m s | , water levels above datum (m) and 
concentrations (kg m~ ) for a specified point in the model 
grid. The input data are presented as a verification of proper 
input data manipulation i this is useful in tracing data 
errors. 



Output produced by routine FRITE is given in Figure A-3. This 
is produced at specified intervals (NFRITE, Table A-2, Card 5) 
and consists of velocities (V and components) in mm s" , 
water levels above datum (mm) and concentrations (mg l" ) 
across the entire grid line by line. For a grid more than 20 
columns wide, more than one print line is produced per grid 
row. Lines containing no non-zero values are deleted. 

Tape output on file 10 consists of groups of five records of 
data representing U and V velocities, water levels above datum 
calculated at half and full time steps and concentrations 
calculated at full time steps. This output is written at 
selected intervals (NTAPEj Table A-2, Card 6) and is used as 
input to cross correlation and plotting routines. 
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TftBIE Al: Files Used by Two-diiiiaisional Numerical Model Program 



File* 
Name 



Number of WDriis 

of Data 
per Itecord 



Logical* 

BBCord 

Laigth 



Block* 
Size 



Description 



I 



FTlOFOOl 

PT20F001 
FT25F001 

FT30F001 
FT35F001 

FT36F001 

FT45F001 

FT46F001 

FT99F001 

'Notes : 



500# 

1 
16 

80 

1 

1## 
1 

1## 
500# 



2004# 

m 



324 
8 

8## 

I 

8## 

20041 



6016# 

1204 
1704 

3244 
804 

eo4## 

804## 
6016# 



Output of velocities, water levels and conc^i- 
trations at selected tiine steps 

Wind ^)eed and direction 

2 
Source flows (dividai by Ax where Ax is the size 

of grid square) 

Source concentrations 

Water ^aeed (signed scalar) at opoi boundary on 
left-tiand side of model grid (row boundary) 

Dissolved- solids concentration corresponding to 
above speed 

Water speed (signed scalar) at open boundary on 
top of model grid (ooluim boundary) 

Dissolved-solids concentration corresponding to 

above speed 

Input of velocities, water levels and ooncentrations 
frcm previous model nm 



*Peculiar to ffl4 360-370 installations. Nuniber wi1±iin file name is valid for any installation, and is 
used in text to refer to that file. 

fPeculiar to Toronto Harbour and must be changed acoording to grid size of area being modelled. 

##For one ocmponeit in concentration model,' must be increased for additional ocnponentB being modelled. 



TART.K j\2: Numerical MDdel Input Card Data 



Card Oolmim 


J Format 
12 


Variable Name 


9-10 


im 


19-20 


12 


NC 


29-30 


12 


IMM 



37-42 



16 



njp' 



67-70 



14 



ngiTIL 






37-42 



16 



N0SC 



37-42 



16 



NSTOP 



DescriptJon 
NuriDer of rcws in model ^id 
Nijiiber of oolums in model grid 
Nuntoer of GOtponients being HiDdelled 



- vAmn model started f rem rest 

= 1 for ipdate nm with one set of input data on file 99 
>1 for i^date run with INITIL/mAPE sets of irput data 
on file 99 (see below) 

Nmtoer of time steps modelled in previous run (not 
used if IUP=0) 



Nimiber of time steps for "spin-ip" of model to initial 
condition (winds are held constant ^Aile source and 
bomda^ flows are interpolated from zero to first 
value m this period) 



Total niitiDer of modelled steps, incliding "spin-ii)" and 
previously modelled steps 



TmiM Mz Oontinued 



Card # card Ooluims Format Variable Wmm 



37-42 



16 



NST^C 



67-70 



14 



nfrue 



37-42 



16 



NTffiE 



67-70 



14 



NSTEPS 



I 



i 



M' 



37-42 
67-70 
31-70 

31-38 
56-65 

21-30 
46-55 
63-70 



16 


DFrr.T 


14 


DET.X 


4F10.0 


MMGN 


Fe.2 


Gmv 


F10.8 


FC0R 


F10.8 


•It 


F10.8 


mm 


F8.2 


KH0 



teecription 

Nuraber of hydro%namic iTDdel time steps for each 
mass balance tiine step. 

Int^val (in time steps) at wMeh output from sub- 
routine FKITE is desired. 

Interval (in time steps) at which model velocities, 
wat^r levels and concentrations are writtai to out- 
put tape (file 10) , 

Interval (in time steps) between successive records 
of real-time input data on disk (files 20-46) . 

Tine step length (in seconds) for hydrodynamic model. 

«^id space step length (in neters) 

Four values of Manning's N (bottom roughness) (Punch 
deciinal point for this and all F-fonmt variables) . 

Acceleration due to gravity (9.815 m s ) 

Coriolis force (0.00010655 s~^ at latitude of Lake 
Ciitario) . 

Wind stress coefficient (dimensionless) . 



Densi^ of air (1.201 kg m ^). 



-3, 



Density of water (999.579 kg m ) . 



TABLE A2: ODntiiiued 





card # 


Card Colinms 


Format 


Variable Name 




U 


36-45 


Flo. 5 


DISPW 




a 


31-35 


15 


NSiC 






65-70 


15 


NI^W 




13 


51-52 


m 


NS0R 




li 


1-64 


1614 


(R0W(I), 00] 




m 


46-61 


F6,0i 


»R 




M 


1-2 


m 


T.TMTT 


1 




17-64 


12(2F2.0) 


CwmDA(i) ,m 



m 



1-3 



4-63 



II 
2013 



Descrip'tion 

2 -1 
Diversion oDef f icient due to wind (m s ) . 

Nijrriber of seoonds between wind readings. 

Nuntter of wind readings (^eed and direction) per 
wind card, 

Niirber of sources and sinks (present rnaxiimiii 16) . 

(R0W(I), q0L(I) , 1=1, NS0R)Rcw and oolumn number of each sour^ and sink (continue 

on another card if needed) . 

Wind direction constant (is angle from true north to 

nodel Y axis measured ooimterclodcwise) . 

Punch "15" (identifies wind card) , 



Wind direction and speed (mi/hr) 

Wind direction according to following 2-digit codes: 
(frOTi) 

11 - N 33 - E 55 - S 77 - W 

12 - NNE 34 - ESE 56 - SSW 78 - Vim 

22 - NE 44 - SE 66 - SW 88 - N 

23 - ENE 45 - SSE! 67 - WSW 89 - tlW 

Rspeat this card for every 12 hr of mnd data required, 
and follo^j with a card which has col. 1-2 blank. 



NCHMIR 



lOffiNR 



Number of open boundary points, on left-hand Mge of 
.raodel 'grid. 

Rcw nureber of cpen boundary points (NCHMR of these) . 



TJfflLE A2: Continiied 



Card # 


Card Coliims 


FO'rmat 


Vari^le Narae 


18 


1-3 


13 


NOMTC 




4-63 


2013 


IjCHMC 


m 


1-5 


15 


NSFCItti 


m 


1-10 


FIO.O 


DT.TM 




21-80 


6F10.0 


■le 



24 



iS 



1-5 



1-10 
21-80 

56-60 



I 
2-70 

56-60 



15 



FIO.O 
6F10.0 

15 



NSiaCi 



DLM 
WLl 



10(F6.0,1X) FI|2ll-a (I) , 1=0, NS0R) 



m 



NSECC 



Descriptioii 

NuiiDer of open boundary points at top 'edge of model 'grid. 

Golum mjitoer of open boiJndary points (NCHM^C of these) , 

Niiriber of seconds between card iiputs of left-hand edge 
opei boundary currents. If these data are specified on an 
input tape, set MSEQ€j = and emit cards 20-22, 

Any nonzero nimber (ideitifies current data) , 

Each 'Card has NCIENR values of left-hand edge open boun- 
dary curraits for one point in time. Repeat every 
teECWL seconds as required, and end with a blank card. 

Nimber of seconds between card inputs of top edge open 
boimdary currents. If these data are specified on an 
input tape, set NSBOtt. - and omit card 22. 

Identical to card #20, open bounda^ current data for 
top edge. 

NurriDer of seoonds beito/een card inputs of source flcM data. 
If these data are specified in an input tape, set NSBCF = 
and emit cards 24 and 25, 

"F" to identify flows. 

1 -1 
flow (m s ) for each source, Ufo cards if NS0R >10. These 

cards are repeated every NSB^" seconds as required and 

temdnated with a blank card. 

Nuifcier of seoonds between card inputs of source concen- 
trations data. If ttese data are specified in an input 
tape, set NSECC = and oniit cards 26. 



TP B JE A2: Oontinued 



Card # Card CoIuehs Pomet Variable Name 



26 



27 



28 



29 



30 



31 



32 



1 
2-70 



Jb-75 



1-50 



1-50 



1-79 



1-50 



1-6 



M MilM 

10(F6.0,1X) (C0NC1(I), 1=1, NS0R) 

7511 m 



a'10.2 



5F10,2 



AMB 



20 (F3. 0, IX) HZ 



5F10.5 



213 



zi^ac 



ID0ISf,JD0I^ 



Description 

"C" to identify ooncentrations. 

Concentrations (ng/l) for eacih souroe, two cards if 
NS0R>1O. These cards are repeated every N^CF seoonds 
as required and temdiiated with a blank card. Hiere 
must be the same nmiiber of cards 24 and 25. 

Land-^water description, NR of these cards, each with 

NC oolumns used. 

4 if land point. 

Blank, 1, 2, 3 for water point, increasing ntmiDer for 

mcreasing value of Manning's N at that grid point 

(aooording to the four N values specified on card #8) . 

MiDient (initial) concentration (rag/l) for each of 
JMM. iiCM3elled ocraponents. 

Lake (outside open boundary) concentrations (mg/l) for 

^di of IMAX modelled conponents, Not used if NSEOCfO 
but card must be present. 

Water depths in feet, IM of these cards, each with NC 
depth values across (me row. If NC>20, each row requires 
2 or more of these cards, according to value of NC, and 
there are NR such sets of cards. 

I^actian matrijc for ip to niax carponents. Zero values 
for conservative siJDstances. 

Rcw and oolunm number of a point for y^ch velocity, 
water level and concentration data mil be printed at 
each tJine step. 



TABIE A3 : Input Data to Heal-Time Cpen Boundaxy Program 
Card # ^Card Columns PO'rmat Variable Naie 





i. 


1-2 


12 




¥!0m 






3=4 


It 




JTHMS 






11-20 


FIO, 





TBC- 






21-30 


FIO, 





TMQ 




1 


1-2 


12 




y0m. 






3-4 


12 




TDME 






5=6 


12 




3^ 


T 


1 


1-3 


13 




li^SC 


to 


' 












4 


1-2 


11 




hUNC 






3-4 


12 




T^mjQ 






5-6 


a 




Mmn? 




t 


1-10 


FIO, 





GmwHii . 






11-20 


FlO. 





TE6MM' 



Descriptiort* 

Start iTDni±i of modelling period (01- Jan, 0'2-Feb, etc,) 

Start day of modelling period* 

Start time for current meter 'data input to model. 

Start time for water 'quality data input to model. 

End month of modelling period,. 
End day of modelling period. 

Year of moielling period,. 

Nlmiber of "spin-i^"' steps,, 

Curreit neter recording interval (min. ) 

Cheriistry 'meter recording interval, (min,, ) 

Model tine length between succeesive inputs of real tiine 
data (min. ) 

Cotpass orientation from which, ,positve (,into irodelled 

area) curT'ente f,low, 

t^axiimm dissolved solids (conductivity *0.65) conoentratlon 
for positive currents If •curr'ent is positive and TDS > 

mamx, TOS is set to 0'.9 * TDSMX. 



TPMJE A4.1 .Input Data tO' ]teal-T,iiie SouT'Ce-Sink 'Progran 



Card # 


Card Coliiwis 


FO'iiTet 


Variable Nane 


1,2,3 










4 


1-10 


FlOi. 


,0 


DKr.X 


5 


1-2 


12 




:iMtTEM 




3-« 


12 




NC^ST 




5-6 


12 




NVAR 


6 


1-2 


12 




MTOC 




3-4 


12 




?tDNQ 


f 


1-2 


12 




MM/ 




3-4 


12 




MINT 



1-10 



F10.0 



CSECT' 



11-20 



F10.0 



CHMBL 



DescriptJQn 

Identical to' cards 'dfis«^iirf in Table A3. 

Grid spaae step length (in netres) . 

Nuitoer of :nietiered souroes (■current .and chaiii,stry meter' 
data avail^le) . 

MiMber of csonstant sources (only average 'data a^ailab'le) . 

mmber of slowly varying sources (have data read in frcm 

'Cards at a specified interval MTOf) 

Current meter reoording interval (min.) , 

'Chemistry meter recordiiig interval (min. ) . This card is 

repeated for each metered soi^*^. 

Nuiber 'Of tours belanjeen readijigs O'f elofing varying sources. 

Itodel tine length between successive .inputs 'Of real time 
'data (min, ) .. 



Cross-section area of water channel , slip or outfall at 
curreit meter location. 

Cofrpass 'Orientation from whi'Ch positive (into modelled area) 
currents flow. ''Itiis card is repeated for 'each metered source, 



TABIiE A4i continued 



Card # 


Card Coluims 


Format 


Variable Nane 


9 


1-72 


12F6.0 


EU2lf 


m 


1-72 


12F6.0 


cm^ 


u 


1-72 


12P6.0 


FL0 



m 



1-72 



12F6.0 



QgSm/ 



Description 

3 
Flow speeds (m /s) of constant sources (Negative for sinks) 

ConcQitrations (mg/l) of ODnstant sources. 

Flow speeds (m /s) of slowly varying sources (negative for 
sinks) . 

GoncQitrations (nng/l) of slowly varYing sources. Qie 
reading on each of cards 11 and 12 every EfflW hours, 
go to entra cards if necessary. JRepeat cards 11 and 
12 for each slowly varying source. If no slowly varying 
source is present, emit these cards. 



FIGURE Al 
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^5. 
























LCR'Ni :i 


300 


UMH'ft 
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10 


'5 




13 


l§ 


22 
2? 


2i 
2? 


















S 


5 


1,2 


ifl 


1« 


7 


17 


20' 


17 


til 12 














7 


13 


13 


16 


31 


2? 


27 


tfe 


2fe 


24 


27 


1 4 9 










^ 


ife 


14 


24 


2 a 


29 


31 


27 


30 


26 


27 


il 


27 


14 










S 


la 


25 


26 


I'f 


2-? 


33 


ift 


2fl 


32 


27 


27 


2§ 










% 


11 


,2B 


30 


2*5' 


a« 


27 


11 


16 


ii 


JS 


■ 17 


27 


20 
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13 


i'^ 


2« 


31 


Si 


31, 


3S 


jp 


50 
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27 


SI 










s 


in 


17 


i'S 


2« 


2B 


2fl 


27 


3'0 


2"5 
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27 


22 








h 


a, 


ifli 


22' 


2'h 
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2« 


27 


ib 


fa 
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27 


27 


20 










S 


Q 


^8 


?a 


^■h 


2 ft 


27 


26 


2h 


24 
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2 7 


27 
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7 


i? 


?3 


?6 


27 


2'5 


27 


ih 
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27 
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% 


I? 


ga 


2'^. 


26 


27 


2S 


ZM 


27 


27 


S'T 
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1 1 


?! 


23 


2'5 


,1'S 


21 


27 


27 


27 


27 
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% 


Q 


21 


21 


21 


2' 3 
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37 
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27 


27 














S 


qi 


31 


20 


Si 
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,27 
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27 


27 
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19 


20 
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27 


27 


27 


27 
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Ife 


19 
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27 


27 


27 


27 


27 
















^. 


|.g 


27 


27 


27 


27 


27 


27 


■2? 


27 






, 






i^ 


p>m 


?.<^ 


?Q 


29 


27 


27 


27 


27 


27 


27 


27 


27 












m 


i^ 


29 


m 


?q' 


2' 7 


27 


27 


27 


27 


27 


27 


27 


















13 


11 


|'7 
20 
20 
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27 
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27 

2fe 
2ft 
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Ifl 


2? 
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FIGURE A2 



NUMERICAL MODEL OUTPUT BY SUB.ROUTINE HARBOR ^ 



T I t4i 


TOTAL 




« INC 


CHANNEL 


F L UW 1 


'SQ^UWCfi 1 


SOURCE 


NSOK 


u 


V 


WATER 


CONO 


" 


31f -> 


VilLUME 




5PECJ 


ID I H 


uum 


'COL 


FLQ'W 


CniNC 


FLOiJi 


CCNC 


13, 6 


IS, 6 


LEVEL 


1 8 » 6 


1 


1 3 (,. 4 \ 


0. 24? 7D 


0'7 


7 . y 2 


62.4 


.00 9 


.048 


J2.30 


'0 . 3 1 


-50,00 


0, 


0.0051 


-0.0 062 


0,46 22 5 


0.23183 


1 J, ^it .^ 


0' 8 242 aD 


7 


7. 92 


&2 .3 


0.00 9 


0,04 7 


32.32 


. J 1 


-50 ,00 


0.0 


- 05 7 


-0.0 049 


0,46 315 


0.23183 


I 


1 -i'H-.J 


.i'4,2QD 


7 


7.92 


62, 1 


O',0 09 


0,047 


32,34 


0,310 


-50. 


0-0 


0.0065 


-0.0 034 


0,46 428 


0.23190 


1 


1 3 4 -1 ^ 


0- 34 3 0'iD 


0'7 


7 . 9 J 


O2.0 


.00 9 


0.0 46 


32.36 


. J 1 


-50-00 


0,0 


0.0071 


-0.OOP4 


0-46 555 


0.23 190 


I 


1 3 4 4 i 


0'. 243 ID 


07 


7*93 


61.9 


0,00 9 


0.046 


.12.36 


0.31 


-50.00 


.0 


0.0075 


-0.0021 


0.46 733 


0-231,90 


-< 


^ 1J446 
1 J44 7 
1344'=? 

1 3 4 4 t:J 


0.243 2D 
. 243i'D 


'0 7 
07 


7.93 

7.') 3 


o 1 , .1 

o 1 . o 


0,009 
.009 


0.046 
0.045 


32.4 
32.42 


0,309 
0,309 


-50. 00 
-50.0 


0.0 
0. 


0.0 077 
0.0076 


-Q«0020 
-0,00 18 


0.4 '6 971 
0-47205 


0» 2,3204 
0.23204 


\ 


#243 30 


C7 


7 . 14 


6 1 , 5 


0. 009 


0.045 


32.44 


0,30 9 


-50. 00 


0,0 


0,0 072 


-0.0021 


0.47390 


0- 23204 




. 24 3 41) 


07 


7.94 


6 1.4 


0.00 9 


0.045 


32.4 5 


0,309 


-50.0 


0. 


0.0 069 


-0.0032 


0.47551 


'0., 232 17 




I 3 '!f 5 


m 243 50 


QJ 


7,94 


'6 1 . 3 


0.009 


0.044 


32,47 


0-308 


-50. 00 


0,0 


0.0070 


-0.0 041 


0.47679 


'D.2321 7 




1 3 4 5 J 


« 243 6D 


01 


7 .9'4 


61.1 


0,009 


,0 44 


32.49 


0,30 8 


-50.0 


0.0 


0.0075 


-0.0048 


0,47741 


0.2 321 7 




i34'-">': 


0. 243 7D 


t>7 


7,95 


61,0 


. 00 9 


.043 


32.51 


.308 


-50. 


0. 


0.0083 


-0.0056 


0-47788 


0, 23222 




13 4 > ^ 


0.24 3 8D 


07 


7.95 


6 0.9 


. D 9 


0.04 3 


32,53 


0,30 8 


-50.00 


0,0 
0,0 


0.0093 

0.0103 


-0-0066 
-0-0074 


0-47 869 
0.47 942 


0. 23222 

0,23,222 




l34'^•+ 


0„243'=>ib 


07 


7 .95 


6 , 8 


0.010 


0.043 


32,55 


0.30 7 


-50. 


I 3 4 o -> 


• 243 9D 


07 


7,95 


6 0,6 


0.010 


,,04 2 


32.57 


0,30 7 


-50. 


0.0 


0.0 109 


-0.0 084 


0,, 4 8 008 


0-23227 




1 3 4 b 6 


, 244 01) 


07 


7 . 96 


60. 5 


, 1, 


0.04 2 


32 .59 


0.307 


-50.0 


. 


0.0 108 


-0-0 096 


0-48125 


0.2 3227 




13 4?; ^ 


,244 lO 


Of 


7.96 


6 . 4 


0.010 


0.042 


3 2,61 


0„307 


-50 . 


0,0. 


0.0103 


-0.0 105 


0-48 24 1 


0.2322 7 




13 4 5 


0. 2442D 


7 


7, 9o 


6 0.3 


U . 1 u 


0.041 


32.62 


0.30 6 


-50. 


0.0 


O.OO'QS 


-0.01 14 


0,48 258 


0.23231 




J 3 4 5 1 


a 2443D 


07 


7.96 


60. 1 


0.01 


0.041 


3 2.64 


0,30 6 


-50,0 


. 


0.008 5 


-0.0127 


0.48 189 


0.23231 




13460 


0, 244 40 


07 


7,97 


60,0 


0.01 


0.040 


32.66 


0.306 


-50.0 


0.0 


0-0076 


-0.0145 


0,48 08 7 


0-23231 




1^461 


0,2 44 40' 
0' • £44 50 


07 
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7,97 
" ■ 7,97 


59,9 
59* 8 


0.010 


0,040 


32.68 


0,306 


-50,00 
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0-0 


0.0073 


-O.0l'$.2 


Of 47 942 


0«23££2 




134 6.^ 


0.010 


0.040 


32.70 


0.30 5 


-50. 00 


0,0 072 


-0,0181 


0.47756 


0,23222 




1 346 1 


0'« 244 60 


07 


7 ,97 


59 . 6 


0.010 


0#039 


32.72 


0.305 


-50,00 


0.0 


0*00?5 


-0,020'5 


0.476 10 


0.23222 




13464 


Om 244 7D 


07 


7.9B 


59.5 


0.010 


0.039 


32.74 


0.305 


-50.00 


0.0 


O'.OOBl 


-0,0228 


0.4 7 559 


. 232 1 '5 




1 3 46 5 


0»244fl.D 


07 


7,98 


59.4 


0.010 


0,0 38 


32.76 


0,305 


-50,00 


0.0 


0*0092 


-0.0,241 
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1 3 4 6 D 


» 2 44 8D 


7 


7.98 


59.3 
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0.038 


32, 7B 


0.305 


-50.00 


0,0 


0.0104 


-0,024« 


0-47563 


0,2321 '5 




1 346 7 


#244 90 


07 


7.98 


59. 1 


0.01 


0,038 


32.79 


0.JQ4 


-50.00 


. 


0.0113 


-0-0 254 


0.47'SS4 


0,23214 
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. 1 0' 
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. 1 1 a 
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0.23214 
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« I 
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-0.0231 


0.4? 41 4 


0-23214 




1 347(5 


" a .res go 


0? 


?«'09 


58. d 


0.010 


0.037 


32.85 


0.30 4 


-50. OO 


. 


. 1 1 e 


-0-021 


0-47 196 


0, 23'2@2 




13471 


0. 245 2'D 


07 


7 #99 


58. 6 
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32-87 
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-50,0 


0.0 


0-0113 


'-0.0193 


0.46996 


0-23202 
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0# 24530' 


7 


8,00 


58,5 


0.010 


0.033 


32.89 


0,303 


-50.00 


0.0 


. 1 '0 7 


-0.0 168 


0,46 866 


0.2 32 02 




13 4 7 3 


. 245 40 


7 


8. '00 


58.4 


0,010 


0,031 


32.91 


0.30 3 


-50.00 


0,0 


0.0106 


-0.0 136 


0, 46 745 


0-23192, 
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• 245 4D 
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S.O'O 
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32.93 


0.303 


-50,0 


. 
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I 34 7'-i 


.245 70 


7 


8.0 1 


57.6 


0.010 


0,019 


33.02 


0.301 


-50.0 


0.0 


0.0097 


-0.0 034 


0,47 056 


0.23201 




i 3 4 (3 


0.245 7D 


7 


8.02 


57.5 


0,011 


« 1 7 


33.04 


0-301 


-5 0.0 


0,0 


0,0082 


-0.0028 


0-47 204 


0.23201 




1 3481 


,245 80 


7 


8,0 2 


57 , 4 


0.011 


0.016 


33.06 


0.301 


-50,00 
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